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ABSTRACT 

We propose that the decline in the near-IR flux from the massive binary sys- 
tem rj Carinae during the spectroscopic event might be explained by accreted 
mass that absorbs the radiation from the secondary star, and by that reduces 
the heating of the dust that is responsible for the near-IR emission. This bi- 
nary system has an orbital period of 2024 days and eccentricity of e ~ 0.9. The 
emission in several bands declines for several weeks near every periastron pas- 
sages, in what is termed the spectroscopic event. In the accretion model for the 
spectroscopic event the secondary star accretes mass from the primary's wind 
for ~ 10 weeks near every periastron passage. The mass is accreted mainly in 
the equatorial plane. The disk and its wind block the secondary's radiation from 
heating dust that does not reside within narrow cones along the symmetry axis. 
This, we propose, might explain the decline in the near-IR flux occurring at the 
beginning of each spectroscopic event. We also argue that the increase in the 
near-IR prior to the event might be accounted for by enhanced hot (T ~ 1700 K) 
dust formation in the collision region of the winds from the two stars. This dust 
resides within ~ 60° from the equatorial plane, and most of it cannot be heated 
by the secondary during the accretion phase. 

Subject headings: (stars:) binaries: general— stars: mass loss— stars: winds, 
outflows— stars: individual (77 Carinae) 



1. INTRODUCTION 

The luminous blue variable (LBV) Eta Carinae has been observed and widely studied 
for almost 170 years, but the existence of its binary companion (Damineli 1996) wasn't 
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fully accepted until about ten years ago. The two winds blown by the massive stellar binary 
system are major players in the 5.54 yrs periodicity, as was recognized in the past particularly 
for the X-ray behavior (Corcoran 2005; Pittard & Corcoran 2002; Akashi et al. 2006). The 
importance of the secondary star and the colliding winds in explaining the behavior at other 
wavelengths has emerged only recently (Abraham et al. 2005; Ishibashi et al. 1999; Damineli 
et al. 2000; Corcoran et al. 2001; 2004; Hillier et al. 2001; Pittard & Corcoran 2002; Smith 
et al. 2004; Steiner & Damineli 2004; Verner et al. 2005; van Genderen et al. 2006; van 
Genderen k Sterken 2007; Soker 2007; Kashi & Soker 2007a, b; Soker & Behar 2006; Behar 
et al. 2007). 

LBVs are massive stars that are believed to be in a rapid and unstable evolutionary 
phase in which tens of solar mass of material are ejected into the interstellar medium over 
a relatively short period of time, possibly in eruptive events (Smith & Owocki 2006; Smith 
2007). Such a twenty year long outburst, known as the 'Great Eruption', was experienced 
by i] Car in 1843-1863 (Davidson & Humphreys 1997), accompanied by ejection of a mass 
of > 12M (Smith et al. 2003a; Smith & Owocki 2006; Smith 2007). This ejecta now 
form a reflecting bipolar nebula— the 'Homunculus'— which has a dense equatorial gas, and 
dense clumps above and below the equatorial plane that obscure the central star at many 
wavelengths. Another outburst, the 'Lesser Eruption' occurred between 1887 and 1895, and 
ejected material mainly in the polar direction, but also in the equatorial direction (Smith 
2005; Smith 2005 & Gehrz 1998; Humphreys et al. 1999). 

The ejected mass has cooled since the Great Eruption, and formed ~ 0.1M© of dust, 
which obscures the central star, reflects light toward us in other directions, and emits in the 
IR band (Smith et al. 2003b). The IR luminosity is L IR ~ 5 x 10 6 L Q (Cox et al. 1995; 
Smith et al. 2003b). According to estimates made by Davidson & Humphreys (1997), this 
IR emission comprises ~ 90% of the luminosity emitted by r\ Car. 

The light curve of r\ Car is characterized by a roughly constant luminosity at all bands, 
although with fluctuations, for most of the 5.54 yrs orbital period (Whitelock et al. 2004; 
Falceta-Goncalves et al. 2005; Abraham et al. 2005; Corcoran 2005; Duncan & White 2003; 
Nielsen et al. 2007). Large variations, typically with slow increase and then rapid decrease 
that follows by a low state that typically lasts tens of weeks, occur every 2024 days in what 
is termed the spectroscopic event, after the disappearance of highly ionized lines (Damineli 
1996; Damineli et al. 2000). The event is thought to occur near periastron passage, where 
the distance between the two stars is r = 1.66 AU (for the binary parameters used by us). 
In the last periastron passage in 2003.5, for example, the sharp X-ray drop occurred in 
JD2452819.6 (Corcoran 2005). 

In this paper we address the behavior of the near-IR light-curve before and during the 
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spectroscopic event, attributing the sharp decline in the near-IR to an accretion event. This 
is the 5 th paper in a series of papers aiming at understanding some of the processes during 
the spectroscopic event by a ~ 10 weeks long accretion event onto the secondary star. The 
accretion event was proposed to explain the X-ray behavior (Soker 2005; Akashi et al. 2006; 
Soker & Behar 2006), and was applied to the visible band as well (Soker & Behar 2006; 
Kashi & Soker 2007b; Soker 2007). 

2. THE INFRARED BEHAVIOR NEAR PERIASTRON PASSAGE 

2.1. The light curve in the near-IR 

In this section we will describe several possible explanations for the near-IR behavior, 
which we however, find to be problematic. We start by describing the behavior of the near-IR 
in the JHKL bands, from Whitelock et al. (2004; Whitelock, P. 2007 private communication) 
who presented a detailed near-IR survey of the JHKL bands, which was taken for a period 
of almost 40 years. The flux of the last cycle (the 2003.5 event) can be seen in Fig. [TJ The 
main characteristics of the last cycle are as follows: 

1. Typical flux values. The longer the wavelength the larger are the average quiescent 
flux (the flux during most of the time when the system is away from periastron) and 
the variations in the flux. The average quiescent flux in J is ~ 160 Jy, while that in L 
is ~ 1440 Jy. 

2. Rapid increase. Approximately ~ 150 days before the event a rapid increase in the 
flux in all wavelengths begins. In J and L the fluxes reach peak values which are 
AFj ~ 40 Jy and AFl ~ 170 Jy above the average quiescent flux. 

3. The event. During the event the hard X-ray emission drops to ~ 1% of the flux before 
minimum and stays close to that level for about 70 days (Corcoran 2005; Hamaguchi 
et al. 2007). The IR event is preceded by a brightening at all near-IR wavelengths, as 
is the case for the X-ray light curve. However, in marked contrast to the X-ray light 
curve the near-IR fluxes decrease during the event by only ~ 15 — 24%. Only in the 
L band the light curve is flat-bottom, and only for ~ 20 days. While the oscillations 
in the X-ray flux make it difficult to pinpoint the time of ingress, it is clearly earlier 
than it is at IR wavelengths. Within the near-IR light curves ingress starts earliest at 
L and latest at J where it is more or less coincident with the well-determined start of 
X-ray flat minimum; the delay between L and J is ~ 7 days. The K minimum occurs 
at X-ray phase 0.011 (~ 22 days after the beginning of the flat X-ray minimum). We 
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will attribute the delay of the J ingress to the formation of extra hot dust close to the 
secondary star near periastron passages. 

4. Typical variation amplitude. Both the relative and absolute decline (the V-shaped 
part of the light curves) during the event are larger for longer wavelengths, being 
AFj ~ —35 Jy, and AFl ~ —400 Jy from peak to minimum. In the L band the 
minimum flux in the event is below the average quiescent flux. 

5. Recover. After the sharp decline, the L band recovers to more or less its average 
quiescent flux, while in the JHK bands the fluxes recover to values above the average 
quiescent flux, and only after the event (lasting ~ 10 weeks) they decline slowly to the 
quiescent level. 



2.2. Problems with a free-free emission model 

Whitelock et al. (2004) argued that the main contribution to the JHKL flux variations 
are due to free-free emission. We now show that this explanation is problematic since free- 
free emission cannot provide enough flux to explain the observations. We assume that the 
recombination rate R = as J n e nidV, is in equilibrium with the ionizing photon emitted by 
the secondary star per unit time $2 — 2.5 — 4.4 x 10 49 s _1 (see Kashi & Soker 2007a). Here 
n e , rii, and dV are the electron number density, ions number density and a unit volume, 
respectively. Using this assumption to substitute for f n e riidV in the expression for the 
free-free emission in the near-IR (Rybicki & Lightmann 2004) we find the contribution of 
free-free emission to the observed flux 
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We took T = 10 4 K in equation ([I]) (for other values of T the emission will be about the 
same or lower), and D n = 2.3 kpc as the distance to r\ Car (Davidson & Humphreys 1997). 
We also used gff ~ 1.2 as the Gaunt factor in the near-IR. The expected flux in the near-IR 
as given by equation (JT]) is plotted in Fig. [2j 

Fig. [2] shows that the fluxes the free-free emission can supply, e.g., Fj = 5.8 Jy and 
Fl = 12 Jy, are much smaller than the typical variations in fluxes as can be seen in Fig. [TJ 
The ionization by the primary star was not taken into account. Adding it to the ionizing 
photon rate of the secondary star would increase the flux by 20% at most. Moreover, these 
results were obtained using 100% production efficiency, and the maximum accepted value for 
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Fig. 1.— Left: The JHKL fluxes for the last cycle taken from Whitelock et al. (2004; White- 
lock, P. 2007 private communication), from top to bottom, respectively. Right: Zooming on 
the 2003.5 event. The left vertical lines in both sides of the figure mark the X-ray decline 
(Corcoran 2005). The right vertical lines mark the end of the assumed 10 week-long accretion 
process (Soker 2005; Akashi et al. 2006) as discussed in section [6] 



the recombination rate R, so we expect the free-free contribution to the near-IR flux to be 
even smaller. We conclude that the free-free emission cannot explain the near-IR behavior. 

We emphasize that our calculated changes in the flux near periastron are attributed to 
ionization by the secondary. The changes, according to the free-free model we are examining 
here, are due to the dense gas in the colliding wind region. This gas must be ionized. The 
primary is much cooler than the secondary, and despite being five times as luminous as the 
secondary has a much lower ionizing flux. Fig. [2] shows that the ionizing radiation cannot 
explain the changes. More problematic to the free-free emission model is the fact that both 
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Fig. 2. — The free-free flux in the near-IR as given by equation [TJ 

the quiescent flux and the variations are larger in the L band than in the K band. The 
free-free emission from the primary wind is expected to behave in the opposite manner as 
in the case of the wind of P Cygni (Lamers et al. 1996). Therefore, although the primary 
wind does contribute to the near-IR emission via free-free emission, it cannot be the main 
near-IR source. Indeed, in P Cygni the near-IR flux is two orders of magnitude below that 
in Eta Car (Lamers et al. 1996). 

2.3. Problems with a free-free absorption 

We now show that the free-free absorption cannot lead to sharp decline in the near-IR 
during the event. The absorption can be explained in principle by the dense conically-shaped 
post-shock primary wind. The decline during the event in J is ~ 15% and in L it is ~ 24% 
(section l2.ip . The free- free optical depth in the L band is ~ 8.4 times larger than that in J. 
Therefore, extinction by ~ 15% in the J band implies almost completely extinction of the L 
band. The similar decline in J and L can be solved in principle if the dense material covers 
all (or most) of the hot dust emitting in the J band, and absorbs ~ 15% of it. The same 
dense gas covers and completely absorbs only ~ 24% of the radiation from the warm dust 
responsible for the L band. This requires a special geometry. 

Another problem with such an absorbing model is that there is no good candidate for 
an absorber: 

1. The dense primary wind close to the primary covers a too small angle, considering that 
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the orbital plane is inclined by ~ 45° to the line of sight, and the IR emitting region 
is expected to be large. 

2. The dense post-shock primary wind material (the conical shell) cannot absorb enough 
of the J-band emission. For details see Appendix A. 

3. The ejection of a dense shell as a general explanation for the spectroscopic event of r] 
Car (Zanella et al. 1984) is ruled out by the X-ray behavior (Akashi et al. 2006). 

We therefore regard any model for the sharp decline in the near-IR that is based on a 
dense gas that obscures the IR emitting regions as unlikely. 



Another model that can account in principle for the sharp decline in the near-IR is 
based on a fast decline in dust formation during the event. Although we do argue for 
near-IR emission from dust, an explanation based on a decrease in the dust formation 
rate for the near-IR decline cannot work for rj Car. The decline in the dust formation 
rate occurs because of a decline in the supply of matter from the binary system and/or 
an enhanced emission that prevents the matter from cooling and forming dust. In 
either cases we would expect that the dust residing closer to the binary system will be 
influenced first, namely, the hotter dust. The hotter dust contributes to the shorter 
IR wavelengths, and would imply that the J-band decline occurs before the Z-band 
decline. This is contrary to observation (Whitelock et al. 2004; see Fig. [I] here). 



Cox et al. (1995) explained the IR light curve of i] Car by a modified blackbody 
emission from dust at two temperatures: A cold dust component at 210 K, and a warm dust 
component at 430 K. We are interested in the hot dust, and for that we will not modify 
the fitting of the cold dust as was done by Smith et al. (2003b; see below), but rather only 
add a third hot component, to better fit the JHKL bands. We will use the same two dust 
components and the modified blackbody function as in Cox et al. (1995) 



2.4. Problems with decreasing dust formation rate 



3. 



FITTING THE LIGHTCURVE 
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where B{T) is the usual Planck function, and m = 1 for amorphous Carbon grains (Hilden- 
brand 1983). 
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In choosing the temperature of the hot dust component we note that Smith et al. 
(2007) used hot dust at T = 1660 — 1750 K to explain the emission from SN2006jc. They 
claimed that although the shape of the red continuum excess can be approximately fit by 
a T ~ 2000 K Planck function, optically thin emitting dust grains will have wavelength- 
dependent emissivity and as long as the grains are not larger than about a = 1 /im, the dust 
will typically have emissivity proportional to with — 1 < f3 < —2. This model is similar 
to the modified blackbody function of Cox et al. (1995). The work of Smith et al. (2007) 
suggests that dust can also reach temperatures of T ~ 1700 K in an environment similar 
to that in the colliding winds of rj Car. Following the above discussion, we will take the 
dust in f] Car to be composed of three components: Hot dust at T = 1700 K, warm dust at 
T = 430 K, and cold dust at T = 210 K. We emphasize that we expect the dust in rj Car to 
have a continuous range of temperatures, and therefore the three dust components used by 
us should be regarded as a representative description of the real situation (a 'toy model'). 
Smith et al. (2003b), for example, fit the cold dust with two temperatures of 200 K and 
140 K. We are not interested in the cold dust in this work, and present the cold dust with 
only one temperature. For that we chose a temperature of 210 K. Our fit cannot be used for 
a quantitative study of the cold gas. In addition, different points in the lightcurve are taken 
from different sources (as summarized in Cox et al. 1995). For these reasons our treatment of 
the cold dust is not accurate, and the work of Smith et al. (2003b) that is superior to our fit 
should be consulted when studying the cold dust. This approximate description is sufficient 
to teach us about the basic physical processes relevant to the near-IR emission. Only the 
two hotter components, the hot and the warm, are relevant to the JHKL wavebands. 

We take the average quiescent flux of JHKL discussed in section [2j add it to the data 
collected by Cox et al. (1995), and create a 3- Temperature modified blackbody fit to the 
entire IR range 

F(X) = b 210J B(210 K) + b 430J B(430 K) + b 1700 5(1700 K). (3) 

The coefficients b\ were obtained by the method of least squares adjustment. The results 
are presented in Fig. [3J 

To better understand the behavior of the JHKL emission during the cycle, and in par- 
ticular near the event we fit the flux from Whitelock (2004) only with the two higher tem- 
perature components, and as a function of time 

F(X,t) = b 430 (05(430 K) +b 170 o(£) -8(1700 K). (4) 

We then converted these coefficients to the luminosity of each dust component by integrating 
over wavelength the modified blackbody flux, and multiplying by AirD; t 2 . The luminosity of 
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Fig. 3. — Three temperature modified blackbody fit to the IR. The points for the JHKL 
wavebands (4 shortest wavebands) are the average quiescent values from Whitelock et al. 
(2004), before the beginning of the event. The rest of the points are taken from different 
sources, summarized in Cox et al. (1995). Each modified blackbody is presented by a thin 
black line, and the overall adjusted flux is presented by a thick blue line. 



the warm and hot dust components as function of time, and in units of their respective 
quiescent average, is given in Fig. HI 



4. DUST PROPERTIES 

Using standard references (e.g., Hildenbrand 1983; Krugel 2003; Whittet 2003) we take 
the following steps. We assume spherical dust grains with radius a = 0.25/im and density 
Pd ~ 1.9 g cm -3 (Hildenbrand 1983). The ratio Q/a, where Q is the coefficient of the 
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Fig. 4. — The luminosity of the hot (1700 K; upper red line) and warm (430 K; lower 
blue line) dust components in units of their respective average quiescent value L q . The left 
panel show an extended period, while the right panel zoom in on the 2003.5 event, together 
with the normalized X-Ray flux (from Corcoran 2005). The X-Ray flat minimum begins at 
t = 2819.6. 

modified blackbody emission (equation [2]) , is obtained from figure 8.1 of Krugel (2003), by 
approximating the wavelengths' range 1/im < A < 100/im 



a 



7.93 



(5) 



For our assumed grain radius a ~ 0.25/im (Hildenbrand 1983; Andersen 2007) we find 
Ao — 2/im (see equation [2]). 

We will now consider a spherical dust grain with a radius a, distance from the star, 
and at equilibrium with the stellar radiation. The grain is heated mainly by the UV radiation 
of the star, which has luminosity L*. At the UV band the absorbtion factor for dust is Q ~ 1. 



The cross section of the spherical dust grain is taken as o 



na 



although it is considered 



a bit larger in some models. The dust will emit the radiation isotropically as the modified 
blackbody (equation [2]). The thermal equilibrium is expressed as 



a 



Ana 2 ■ tic / B(T)\- 2 d\. 



(6) 



where T is the dust temperature. After integrating over the modified blackbody, we get 
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where we scaled the stellar luminosity with the assumed secondary's luminosity. The reason 
is as follows. We assume that most of the variation in the hot dust luminosity is from hot 
dust formed in the shocked primary wind as it collides with the secondary wind. The dense 
gas of the primary wind close to the primary, and in the conical shock region formed by 
the wind collision absorb most of the primary radiation (Kashi & Soker 2007a). Cold dust 
which is distributed at larger distances is heated mainly by the primary stellar radiation, and 
possibly by the secondary radiation that was not absorbed by the hot dust. Hence we take 
L* = L 1 + L 2 = 5 x 1O 6 L (Cox et al. 1995; Davidson & Humphreys 2000) for T = 430 K 
and T = 210 K. We take two cases for the hot dust at T = 1700 K : L* = L 2 = 9 x 1O 5 L 
(Verner et al. 2005) as a lower limit, and = Li + L 2 as an upper limit. The lower limit 
gives r d = 26.8 AU, as seen from equation [71 and the upper limit gives r d = 63.1 AU. Van 
Genderen et al. (1994) found r d = 63 AU for T = 1700 K using regular blackbody function 
with totally different parameters which were acceptable at the time. 

We can also evaluate the total mass of the dust. By taking the luminosity that N d (T) 
grains would produce equal to the observed luminosity obtained by integrating the adjusted 
flux of each component (equation [3]) over wavelength 

POO 

L d (T) = c bi5(T)A- 2 dA, (8) 
Jo 

we can find the number of grains from the equation 

poo 

L d (T) = N d (T)47ra 2 ■ tcc / 5(T)A~ 2 d\, (9) 

Jo 

and then find the mass of the dust from 

4-7T 

M d (T) = N d (T)—a 3 Pd . (10) 

Taking p d ~ 1.9 g cm -3 (Hildenbrand 1983) for the grain density, and the following observed 
quiescent luminosity for each dust temperature Ld(1700 K) = 1.1 • 1O 5 L , L^(430 K) = 
7.7- 1O 5 L , and L d (210 K) = 3.9 • 1O 6 L , we find the approximate dust mass (only dust, not 
the inferred total gas mass) M d (1700 K) ~ 7.6 x 1O~ 9 M , M d (430 K) ~ 1.4 x 1O" 5 M , and 
M d (210 K) ~ 8.4 x 1O" 5 M . 

In an attempt to estimate the mass of dust that resides at a distance r d (T) from the 
binary system we assume that gas flows from the binary system at the terminal velocity 
of the primary wind v\ ~ 500km s -1 . The time it would take for the gas to arrive at 
distance r d from the star is tfl ow (^) = r d(T)/vi. For the primary mass loss rate we take 
Mi = 3 x 1O _4 M , and the conical shock to cover ~ 1/4 of the space around the star (see 
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Kashi & Soker 2007a). Assuming the gas to dust ratio is ~ 100, the estimated dust mass 
from flow consideration is 

Mr(T) * « * = 7 x 10- 7 i ;ff U3xl0 4 syr . 1 M e (11) 

In this formula the flow time ta ow (T) depends on the distance of the dust, which in turn 
depends on the dust temperature. We summarize our findings in Table 1. 

We emphasize that the last two columns of Table 1 are calculated for the (wrong) 
assumption that the dust is supplied by the present wind. Clearly, as we explain below 
and is well known (Smith et al. 2003b), the warm and cold dust cannot be supplied by the 
present wind. Therefore, for example, the flow time to the warm dust of ~ 18.6 years is not 
the real flow time. It is the flow time had the dust been supplied by the present wind. The 
hot dust can be supplied by the present wind. Although our "toy model" does not intend to 
give exact results to the distances and ages of the dust, but only to give a general picture, we 
bring and discuss Table 1 here because it emphasizes the differences between the hot dust 
on one hand, and the warm and cold dust components on the other hand. 

From Table 1 we learn that the cold dust, which is responsible for most of the observed 
IR radiation from rj Car (Cox et al. 1995), resides at very large distances from the star, 
compatible with its ejection during the Great Eruption that ended ~ 150 years ago (Smith 
et al. 2003b; Smith & Gehrz 1998). A small part of this dust was ejected during the Lesser 
Eruption ~ 115 years ago, which created circumstellar gas at distances of a few hundred to 
a few thousand AU from the star, mainly in the polar directions, but some in the equatorial 
plane as well (Smith et al. 2005). Indeed, the present regular flow from the primary cannot 
supply the mass required to explain the emission by the cold dust, as M| ow <C M d . As 
stated earlier, our analysis of the cold dust is not accurate, but meant only to present the 
basic differences between the hot and cold dust, e.g., that the cold dust cannot be supplied 
by the present mass loss from r\ Car. A proper analysis of the cold dust, that includes also 
colder dust at 140 K, yields even larger dust mass of ~ O.1M (total gas and dust mass of 
~ 1OM ) in the cold phase (Smith et al. 2003b). The required cold and very cold dust mass 
of ~ O.IMq (a gas mass of ~ 1OM ) could be easily supplied by the mass lost during the 
great eruption (Smith et al. 2003b; Smith 2006; Smith & Ferland 2007). The total luminosity 
of all three components L tot = 4.75 x 10 6 L© is smaller than the assumed luminosity of r\ Car 
L v = 5 x 1O 6 L (Cox et al. 1995). This and our finding that the dust mass can be much 
smaller than that ejected during the great eruption, implies that some radiation escapes 
without being absorbed by dust. Probably, mainly along the polar direction (Smith et al. 
2003a). 

Smith et al. (2003b) have also mapped the dust of rj Carinae. They used three different 



13 



A K) 


r d {AU) 


L d (L Q ) 


M d (M Q ) 


Uow( yr) 


M| OW (M ) 


1700 (lower limit) 


26.8 


1.1 x 10 5 


7.6 x 10~ 9 


0.25 


~ 2 x 10~ 7 


1700 (upper limit) 


63.1 


1.1 x 10 5 


7.6 x 10~ 9 


0.60 


~ 5 x 10~ 7 


430 


1,960 


7.7 x 10 5 


5.3 x 10~ 5 


18.6 


~ 1.5 x 10~ 5 


210 


11,756 


3.9 x 10 6 


9.6 x 10~ 3 


111.5 


~ 8 x 10~ 5 



Table 1: The properties of the three dust components. T( K) is the dust temperature 
assumed in our three dust components model. r d is the distance of the dust to the star, 
calculated by assuming that the hot dust is heated either by the secondary star or both 
stars, while the warm and cold dust components are heated by the primary and some of the 
secondary radiation (primary and secondary luminosities are given in the text). L d is the 
luminosity of each component as calculated from the observed IR flux (equation [8]). M d is 
the dust mass calculated from the observed luminosity by equations © and (|T0|) . The last 
two columns examine whether the present primary wind can supply the dust. Therefore we 
take the outflow velocity to be V\ ~ 500km s -1 . t^ ow = r d /v\ is the flow time of the gas from 
the binary system to the dust location r d . M| ow is the dust mass residing near r d crudely 
calculated from the flow time according to equation pip . The results for the 210 K and 
430 K dust show that the present wind cannot supply the dust. Those components resides 
in the Homoncolus and came from the 19th century eruptions. The 1700 K dust forms close 
to the binary system and can be supplied by the present wind. 
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temperatures (140 K, 200 K, 400 K), in their goal to find the total mass and IR luminosity of 
the Homunculus. For a component at T = 200 K they got luminosity L 20 o = 1.2 x 1O 6 L and 
dust mass M 2 oo = 0.015M o , for a component at T = 400 K they got luminosity L 400 = 5.5 x 
1O 5 L and dust mass M 40 o = 2 x 10~ 4 M Q . Their total IR luminosity was L IR = 4.3 x 1O 6 L 
(Smith et al. 2003b). Although using some different parameters our results are quite close 
to their results. 

From Table 1 we also learn that a large part of the warm dust can be formed in the 
outflow from the primary star, and all the hot dust forms from the primary outflow relatively 
close to the binary star. In the next section we propose that the hot dust component and 
part of the warm dust component are formed in the conically shaped post-shock primary 
wind region. From Table 1 we learn also that if for some reason the outflow from the 
binary system changes its properties, the change in the emission of the hot dust will lag by 
~ 0.25 year = 3 months. This time is longer than the ~ 10 weeks duration of the event. 
This suggests that the changes in the IR flux during the event is not due to changes in the 
outflow that forms the dust. The much slower increase in the IR flux starting ~ 150 days 
before the event can be accounted for by changes in the flow properties, as we explain in the 
next section. We attribute the rapid changes in the IR fluxes to obscuration of the secondary 
star by the mass accreted onto it during the event (sections 6.1 and 7). 



5. HOT DUST FORMATION 

Hot dust is present at all times as is inferred from the J and H bands. It is probably 
formed in the primary wind, both the shocked and unshocked segments. Here we don't con- 
sider the dust formation in the primary wind in quiescence. We only examine the possibility 
that the increase in the near infrared luminosity close to periastron passage is caused by 
enhanced dust formation in the shocked primary wind. 

The two stellar winds collide and go through two respective shock waves, with a contact 
discontinuity between the two post-shock wind regions, as drawn schematically in Fig. [5J 
The post shock primary wind cools fast as it outflows (Kashi & Soker 2007a; Soker 2003; 
Pittard & Corcoran 2002) and forms a thin shell in a conical type shape, which we refer 
to as the conical shell. We expect dust to form in this cool high density conical shell. The 
shocked secondary wind has a very low density and does not cool. No dust will be formed 
in that region. The density in the conical shell depends on several parameters, e.g., the 
properties of the secondary winds (Kashi & Soker 2007a). More influential and relevant to 
the present study is the dependence of density on the primary mass loss rate Mi, on the 
orbital separation r, on the distance from the primary r 1; and on the magnetic field in the 
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primary wind via its influence on the post-shock density (Kashi & Soker 2007a). Typical 
values of the proton number density in the conical shell are n p ~ 10 9 — 10 14 cm -3 (Kashi 
& Soker 2007a). It is interesting to note that Smith et al. (2007) require a density of 
n p ~ 10 10 cm~ 3 to exist in the region where the hot (~ 1600 K) dust is formed in SN 2006jc. 

We examine the ionization structure of the conical shell as a result of the secondary 
radiation. Our motivation is to examine whether extra hot dust can be formed in the conical 
shell as the system approaches periastron, in addition to dust expected to be formed in the 
primary wind. We neglect the ionizing radiation from the primary because most of it is 
assumed to ionize the primary wind close to its origin (Kashi & Soker 2007a). We identify 
three stages of ionization along the orbit: (i) When the orbital separation is large, as is the 
case for most of the orbital period, the density in the conical shell is low, and reaches only 
n ~ 10 9 — 10 10 cm -3 in the densest parts. Therefore, the hot secondary ionizing radiation 
manages to keep the conical shell totally ionized. In that stage dust can be formed only at 
very large distances. Not much hot dust will be formed close to the secondary. The hot dust 
might be formed in the primary wind and be heated by the primary and secondary radiation 
which penetrate the conical shell, (ii) In a calculation described in Appendix B below, we 
find that when the orbital separation decreases to r < 14.3 AU at ~ 170 days before (after) 
periastron, the density of the gas in the conical shell increases such that the cone starts 
to be optically thick to the secondary's UV radiation and parts of the cone begin to be 
neutral. First at large distances, and then closer and closer to the stagnation point along 
the line joining the two stars. (Hi) As the secondary approaches periastron, r < 3.2 AU 
at ~ 20 days before (after) periastron the ionization front is within the conical shell in all 
directions, and the region near the stagnation point becomes neutral, as shown in Figs. [9] 
and [TU1 in Appendix B. 

According to the parameters used here, stage (ii) starts ~ 170 days before periastron 
passage. We propose that the favorable conditions for dust formation lead indeed to dust 
formation close to the secondary star. This hot dust is the explanation for the rapid increase 
in the near-IR starting ~ 150 days before the 2003.5 event. Considering the many uncer- 
tainties, e.g., the value of the magnetic field in the primary wind (see Kashi & Soker 2007a) 
and the many fluctuations in the near-IR light curve (Whitelock et al. 2004), we consider 
satisfactory that we find the favorable conditions for hot dust formation, hence the rapid 
increase in the near-IR flux according to our model, to start ~ 170 days before periastron, 
while the observations show it to start ~ 150 days before periastron passage in the 2003.5 
event. The exact time when the conditions become favorable for dust formation is very 
sensitive to several of the winds parameters, and therefore we cannot predict the exact time 
when the favorable conditions for dust formation will start in future events. Indeed, in the 
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Fig. 5. — Schematic drawing of the collision region of the two stellar winds and the definition 
of several quantities used in the paper. There is an axial symmetry around the line through 
the two stars. The two thick lines represent winds' stream lines. The two shock waves are 
drawn only in the lower half. The post-shock regions of the two winds are hatched. The 
shocked primary wind region is referred to as the 'conical shell'. The dashed line shows the 
accretion column which exists, according to the proposed model, only for ~ 70 — 80 days 
during the accretion period which corresponds to the X-ray minimum and the spectroscopic 
event. 

1992.5 and 1998 events there is no rapid increase in the near-IR flux before the event, but 
rather a gradual increase that seems to be part of the secular variation (Whitelock et al. 
2004). In the 1981 event there is a rapid rise before periastron, while in the 1987 event the 
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rapid rise occurs earlier on, starting ~ 300 days before the event and finishing after several 
tens of days. The secular (gradual) increase in the near-IR flux is observed to occur over 
the last 4 decades (Whitelock et al. 2004; Whitelock, P. 2007, private communication). In 
addition to the gradual increase there are slow large flux fluctuations during the orbital pe- 
riod. These fluctuations and the gradual increase will be treated in a separate paper, where 
the secular increase in the near-IR emission will be attributed to the recovery of both the 
primary (from its high mass loss rate) and the secondary (from its high mass accretion rate) 
from the 'Great Eruption'. 

The orbital separation ~ 170 days before periastron passage is r ~ 14.3 AU. The extra 
hot dust starts to be formed at a distance of r 2 ~ 13.4 AU from the secondary. This is 
compatible with the distance of the hot dust from the secondary r 2 ~ 25 — 60 AU as we 
found in section 0] (see Table 1). 

The hydrogen and helium visible and infrared lines are no problem in explaining near-IR 
emission from dust, since they are not generated in the same place where the dust resides. 
The Hel lines are formed in the acceleration zone of the secondary wind (Kashi and Soker 
2007b) while the dust emitting the near-IR is formed in the outer regions of the conical 
shock. 

Walsh & Ageorges (2000) discuss dust features in rj Car, mainly in the Homunculus. 
The IR spectrum of rj Car is characterized by a peak around 10/im (Mitchell & Robinson, 
1978), indicative of silicate grains (Whittet 2003). Walsh & Ageorges (2000) claim that 
few percents of the visual and JHK radiation is polarized. Polarized radiation is a strong 
indication for non-spherical grains. Magnetic fields make the grains align more or less with 
their long axes parallel to each other, causing polarized radiation scattering (Whittet 2003). 

The existence of non-spherical silicate grains instead of spherical amorphous-carbon 
grains, may slightly change Q(\), because the extinction curve for grains from different 
materials and shape is different (Krugel 2003). Consequently the curve of the modified- 
blackbody function we used (equations [2] and [5]) is changed. These small changes, however, 
cannot cause a serious change in our results. 

6. THE ACCRETION PHASE 

In total, the accretion process lasts ~ 10 weeks (Soker 2005; Akashi et al. 2006). We 
consider two modes of the accretion process; we expect the real accretion process to include 
both. The discussion below is based on previous papers of the series, in particular Soker 
(2005b) and Akashi et al. (2006). 
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6.1. 



General Considerations 



As was explained in previous papers in the series (Soker 2005b; Akashi et al. 2006), 
there are strong arguments supporting an accretion phase lasting ~ 10 weeks near periastron 
passage. The disk has time to establish itself because its Keplerian time tad is much shorter 
than the length of the accretion phase. For a secondary mass of M 2 = 3OM the Keplerian 
(dynamical) time at radius r^d in the disk is tdd — l-4(r^/6O-R ) 3 / 2 week. We took the 
accretion disk outer radius to be about three times its inner radius (the secondary radius). 
The parts of the disk closer to the secondary surface establish themselves in less than a week. 

However, the viscous time scale is longer, and the disk does not reach a steady state. The 
viscosity interaction time to spread material in geometrically thin disks is t v d ~ O.lr 2 ^/ (addC s H) 
(Frank et al. 1985), where H ~ (C s /vK)fdd is the vertical size of the disk, vk is the Keplerian 
velocity at rjd, C s is the sound speed, and add ~ 0.01 — 1 is the viscosity parameter. This 
gives that the viscous time scale is 



For H/r^d — 0.3 and a^d — 0.1, we find the viscous time scale to be somewhat longer than 
the dynamical (Keplerian) time scale. The viscous time scale might become as long as a few 
weeks at the outer boundary of the accretion disk (although the time is much shorter close 
to the secondary surface). In our model, therefore, a disk is formed, but it does not reach 
a complete equilibrium, and does not settle into a fully thin disk. The disk probably has a 
wind as part of the evolution toward equilibrium. The disk and its wind block the radiation 
in most directions, allowing most of the secondary radiation to escape only in narrow cones 
in the polar directions. This is drawn schematically in Figure [6j 

For our purpose any dust residing outside the illuminating cones, i.e., away from the 
polar directions, is practically close to the equatorial plane (even if geometrically speaking it 
is not). In particular, this holds for dust that is/ was formed in the conical shell. Therefore, 
most of the excess hot dust that is being formed in the conical shell close to periastron 
passage will not be heated by the secondary radiation during the accretion phase. The same 
holds for dust that is being formed in the accretion column during the accretion phase. 

The secondary luminosity constitutes ~ 18% of the total luminosity in our model (see 
section H]). Because the secondary is closer to the hot dust and a substantial fraction of 
the primary radiation is blocked from heating the excess hot dust (that was formed in the 
conical shell) by the primary's own dense wind, the secondary share in heating the excess 
hot dust is much larger than its ~ 18% contribution to the bolometric luminosity of the 
system. Therefore, during the accretion phase the heating of the hot dust can potentially be 
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Fig. 6. — Schematic geometry of the hot dust and secondary illumination of the hot dust. 
Some components are drawn only on one side of the equatorial plane. Upper panel: When 
the conical shell, formed by the winds' collision, exists, then the secondary illuminates it in 
all directions. Lower panel: During the accretion phase. The accretion disk and its wind 
block the secondary radiation in most directions; radiation escapes in narrow cones along the 
polar directions. The secondary wind does not exist, or exists only along the polar directions. 
Hot dust is formed in the accretion column. 

reduced more that ~ 20%. As we will discuss below, this is indeed the case. 

6.2. The collapse of the conical shell 

In the first mode of accretion the conical shell close to the stagnation region collapses 
first onto the secondary. The optical depth of the conical shell at three orbital phases 
and as function of direction from the secondary is shown in Fig. [3 For the half space of 
a > 90° close to periastron passage {6 = 0) we take an average value of r — 0.4. As this 
region collapses toward the secondary the optical depth increases as r > Tor2 S /r2, where 
r 2s — 0.5 AU and r 2 are the distances of the shell from the secondary before the collapse 
starts and after, respectively. The reason the opacity increases faster than l/r 2 is that the 
density increases as l/r| as the conical shell collapses, and with it the opacity. We find that 
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t ~ 1 at r 2 ~ 0.3 AU ~ 3R 2 . Because of the non-negligible angular momentum, the mass 
will be accreted mostly from near the equatorial plane (Akashi et al. 2006). 

The stagnation point is at a distance of ~ 0.5 AU from the secondary, but regions away 
from the stagnation point are at larger distances. Taking the last region to collapse to fall 
from ~ 1 AU, we find the free fall time to be ~ 12 days, which is about half the duration of 
the fast decline phase. 
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Fig. 7. — The optical depth of the conical shell as a function of direction a and for different 
orbital angles 9: 9 = 90° and 9 = 60° corresponding to ~ 20 and ~ 9 days before periastron 
passage, and 9 = at periastron passage, a = 180° corresponds to the stagnation point. 



6.3. Accretion from the acceleration zone of the primary wind 

In the second mode we assume that near periastron passage most of the mass lost by the 
primary is accreted by the secondary M acC 2 ~ Mi. We note that the accretion rate can be 
even larger as the primary almost fills its Roche lobe. The primary spin is not synchronized 
with the orbital motion, and it is probably slower, the true primary's Roche lobe is smaller 
than calculated for synchronized systems. For an effective temperature of T c gi = 20, 000 K 
and luminosity of L\ = 4.5 x 1O 6 L we find R\ ~ 180i? Q . Assuming that the primary wind 
velocity v\ = 500 km s" 1 is the escape velocity from the primary star, and taking its mass to 
be Mi = 12OM , we find the same radius Ri ~ 18Oi? . The orbital separation between the 
two stars near periastron is in the range a p ~ 250 — 350i? Q for an eccentricity of 0.93 — 0.9, 
respectively. Namely, near periastron Ri/a p ~ 0.5 — 0.75. The secondary is within the 
acceleration zone of the primary wind, and the accretion rate can be very high. 



Assuming that the mass falls on the secondary at the free fall speed Vff — (2GM2/r2) 1//2 , 
and taking for the secondary mass M 2 = 3OM , the density of the accreted matter is 
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For a temperature of 10 4 K we find the opacity to be k ~ 0.3. Therefore, the optical depth 
is 



The angular momentum of this material is large, and most of it will be accreted from the 
equatorial plane. We have no way to estimate the total duration of the high accretion rate 
phase, and have to assume that it lasts ~ 5 weeks, as the duration of the low IR phase. 

Although the conical shell does not exist during the accretion phase, dust continues to 
form in the accretion column— the dense region formed behind the accreting secondary (see 
Fig. [5] and previous papers in the series). The high density and high opacity near periastron 
passage and during the initial accretion phase implies that dust will be formed close to the 
secondary star. This will be hot dust. Therefore, there are two competing effects at the 
beginning of the accretion phase. On the one hand obscuration of the secondary radiation in 
the equatorial plane reduces dust heating. This causes the early decline in the L band. On 
the other hand dust is formed close to the secondary star. This dust is hot, and contributes 
mainly to the J band. Initially this effect dominates over the obscuration of radiation, and 
therefore the flux in J starts to decline later. The delay between the warm and hot dust is 
seen also in Fig. m where the X-ray light curve is also presented. The obscuration becomes 
the dominant effect when the accretion phase is well developed, when the X-ray drops as 
well. 

The accretion phase continues after the secondary has left the acceleration zone of the 
primary wind, and lasts for a total of ~ 10 weeks. The accreted gas at the second half of 
the accretion phase has angular momentum that is not enough to form an accretion disk. 
However, it is high enough to channel accretion to the equatorial plane; we continue to 
refer to this equatorial accretion as 'disk'. The accretion of gas mostly from the vicinity of 
the equatorial plane during the entire accretion phase results in a collimated polar outflow 
(wind) that after being shocked emits in the hard X-ray, and accounts for the residual X-ray 
emission during the event (Akashi et al. 2006 section 4.3.1). The speed of this polar outflow 
is similar to the secondary wind speed, hence will lead to similar X-ray temperature, as 
found by Hamaguchi et al. (2007; for that the criticism of Hamaguchi et al. of the accretion 
model is not valid). We emphasize that we are referring here to a collimated flow blown by 
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the secondary star. The primary star continues to blow its much more massive wind, which 
is spherical during the event (Smith et al. 2003a). 

This polar outflow will make the optical depth in the polar directions much lower than in 
the equatorial plane. Hence, a large fraction of the secondary radiation will be emitted along 
the polar direction during the accretion phase. Because all the excess hot dust (that is formed 
in the conical shell close to periastron passage) resides within < 60° from the equatorial plane, 
and a large fraction of the hot and warm dust that is formed in the unshocked primary wind 
resides away from the polar directions as well [Smith et al. 2000; cooler dust at T < 200 K 
resides in the polar lobes and larger distances (Smith et al. 2003b; Gomez et al. 2006)], we 
expect the dust emission to decrease as the disk blocks the secondary radiation; this is the 
subject of section As the spectroscopic event proceeds, the secondary's distance from the 
primary increases, and the accretion rate decreases. The polar opening angle is expected 
to increase, and more of the secondary radiation will reach regions closer and closer to the 
equatorial plane, where hot dust resides. 

7. DISCUSSION AND SUMMARY: A MODEL FOR THE NEAR-IR 

LIGHTCURVE 

The main goal of this paper is to explain the fast decrease and then increase in the 
near-IR emission that occurs during the spectroscopic event (Whitelock et al. 2004). For 
that, we now incorporate all the ingredients discussed in previous sections, and try to explain 
the near-IR light curve near periastron passage. 

In section 2 we argued that free-free emission cannot account for the near-IR behavior. 
Basically, the free-free emission cannot supply the required luminosity. We also showed 
that absorption of the near-IR emission cannot account for the decline. In sections [3] and H] 
we explored the dust properties, by building a 'toy model' where the dust is composed of 
three-temperature components. We found there that the hot (1700 K) dust can be supplied 
directly by the present mass loss rate of the primary star. The warm and cold dust, on the 
other hand, must come mainly from the eruptions of the 19th century and even from earlier 
times (Smith et al. 2003b; Gomez et al. 2006). 

In section |5] we argued that as the two stars approach periastron conditions become 
favorable for dust formation in the conical shell— the post shock region of the primary wind 
(see Fig. [5]); This is presented in Figs. [9] and [10]). We attribute the rapid increase in the 
near IR flux starting about five months before the 2003.5 event to enhanced dust formation 
in the conical shell. The conditions are very sensitive to several of the wind parameters of 
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both stars. Indeed, in some past events there was no rapid rise, while in others the rapid rise 
occurred at different orbital phases before periastron (see section [5]). We therefore cannot 
predict the exact time for the favorable dust formation conditions, hence the rapid increase 
in the near-IR flux, to occur in future events. In addition to this sensitivity there is the 
gradual change in the IR flux observed for almost 40 years now, which we attribute in a 
future paper to the recovery of both the primary and secondary from the Great Eruption of 
the 19th century. 

We turn to describe the decline during the spectroscopic event. In previous papers in the 
series the absorption of the secondary ionizing radiation was considered (Soker 2007; Kashi 
& Soker 2007a). Here we are also interested in the near UV and visible bands which are 
responsible for heating the dust. Our findings in section [H] suggest that the mass accretion 
rate near periastron passage is such that the secondary radiation capable of heating the 
dust is almost completely blocked in directions away from the polar directions. Most of the 
radiation will escape in narrow cones in the polar directions, and some will be converted to 
longer wavelengths by the optically thick accreted gas. These longer wavelengths are less 
efficient in heating the dust. We expect some of the hot dust to be formed in the conical shock 
(section EJ), hence to reside within ~ 60° of the equatorial plane. The polar radiation will 
not heat this dust, and the dust will cool down and the emission in the JHKL IR bands will 
drop. Dust at very large distances will be heated by the primary's radiation and secondary's 
polar radiation, so we expect no decline in the far IR. (e.g., in ~ 30/zm). 

Most of the excess hot dust, i.e., the dust that is formed within several months from 
periastron passage, is heated by the secondary star, that accounts for ~ 18% of the bolometric 
luminosity. Therefore, the blockage of the radiation of the secondary will erase almost 
completely the increase in the near IR radiation prior to the decline. Hot dust that is 
formed close to the secondary will partially compensate for this decline. More than that, the 
blocked radiation from the secondary (during the accretion phase) does not heat the dust 
that is formed in the unshocked primary wind. Therefore, the total decline can bring the 
radiation to be below its quiescent level. 

The high density and high opacity near periastron passage and during the initial ac- 
cretion phase implies that dust will be formed close to the secondary star. This will be hot 
dust. Therefore, there are two competing effects at the beginning of the accretion phase. On 
the one hand obscuration of the radiation of the secondary in the equatorial plane reduces 
dust heating. This causes the early decline in the L band. On the other hand dust is formed 
close to the secondary star. This dust is hot, and contributes mainly to the J band. Initially 
this effect dominates over the obscuration of radiation, and therefore the flux in J starts to 
decline later. The delay between the warm and hot dust is seen also in Fig. HI where the 
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X-ray light curve is also presented. The obscuration becomes the dominant effect when the 
accretion phase is well developed, at that stage the X-rays sharply drop as well. 

We attribute the sharp decline in the near-IR emission near periastron passage (the V- 
shaped region in the light curve) to the high mass accretion phase occurring at the beginning 
of the accretion phase (lasting ~ 5 weeks). The Bondi-Hoyle accretion rate of the secondary 
from the primary wind outside its acceleration zone is M acc ~ 1O _6 M (Akashi et al. 2006), 
and the optical depth is very low. Therefore, in the second part of the accretion phase 
(lasting ~ 5 weeks,) when the orbital separation increases, the radiation of the secondary 
manages to heat more and more dust closer to the equatorial plane, mainly hot dust, and 
the near-IR emission recovers. 

The decrease in the heating radiation from the secondary toward low latitudes starts 
when matter starts to be accreted on the secondary. This occurs a few weeks before the event. 
This explains the early small decrease in the L band. At the same time, more dust is formed 
close to the secondary star. This is a hot dust, that more than compensate which more 
than compensates for the decrease in the radiation, hence the flux in the J band continues 
to increase. Eventually, as the accretion phase is developed, there is a rapid decline in all 
bands. 

We suspect that after the high mass accretion rate phase starts, it is mainly the nearby 
dust that formed in the accretion column that is being heated. This is a relatively hot dust. 
The secondary does not heat dust at large distances in the equatorial plane, which emits 
much more in the L band. Therefore, J recovers relatively more than L (in absolute flux L 
changes more). 

At all times the primary continues almost unaffected its role in heating the dust, there- 
fore the decline in the IR is only ~ 10 — 24%. 

We are grateful to Patricia whitelock, Freddy Marang and Francois van Wyk for sending 
us their data plotted in Fig. [TJ We thank Nathan Smith for very helpful comments. This 
research was supported by a grant from the Asher Space Research Institute at the Technion. 
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A. Appendix: Free-Free Absorption Calculations: 

We now show that the dense post-shock primary wind material (the conical shell) cannot 
absorb enough of the J-band emission (for more detail on the absorption by the conical shell 
see Kashi & Soker 2007a). The free- free absorption coefficient at frequency v is given by 
equation (5.18) of Rybicki & Lightmann (2004) 

a f J = 3.7 x 10 8 T-^Z 2 n e niU- 3 

where, as discussed in the section 2.2, T = 10 4 K and gjf ~ 1.2. We calculate the free-free 
optical depth for the J-band as a function of a (the angle between a direction opposite 
to the direction of the primary to a point on the contact discontinuity measured from the 
secondary; see Fig. [TU] below) and 9 (the orbital angle; see Fig. [9] below). Defining d as the 
post shock primary wind depth (see Kashi & Soker 2007a, and Fig. [10] below) we find that 

ri f = 2.2 x 10" 37 n e («, #K(a, 6)d(a, 9), (A2) 

where n e , rii and d are in c.g.s. units. Figure [8] shows The free-free optical depth in the J 
waveband as a function of a and 9. The maximum value of Ty is obtained for a direction 
through the primary (a = 180°) at periastron (9 = 0°). For the point on the conical 
shell along this direction the electron and ion densities are n e ~ 1.9 x 10 12 cm~ 3 , and 
Hi ~ 1.6 x 10 12 cm -3 , respectively, and the geometrical width of the conical shell along this 
direction is d ~ 2.5 x 10 11 cm. The maximum value obtained is therefore ~ 0.17. We 
note that this is an hypothetical maximum, because according to our model the colliding 
winds cone does not exist close to periastron passage. We conclude that the dense post- 
shock primary wind (the conical shell) is optically too thin to absorb enough of the J-band 
emission as seen in observations. 
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Fig. 8. — The free-free optical depth in the J waveband as a function of a and 9 as given by 
equation IA21 a is the angle of a ray through the conical shell measured from the secondary 
(a = 180° toward the primary; see Fig. [10] below), and 9 is the orbital angle (9 = 180° at 
apastron and 9 = 0° at periastron). 



B. Appendix: Ionization and Dust Formation in the Colliding Winds Cone: 

In finding the ionization front location we use our colliding winds model (Kashi & Soker 
2007a) and assume that dust forms in the post-shock primary wind region. We define the 
direction angle a as the angle from a direction opposite to the primary direction to a point 
on the contact discontinuity as measured form the secondary (see Fig. [10]). Segments of 
the conical shell start to be neutral when the total recombination rate per steradian along a 
direction a measured from the secondary becomes equal or larger than the rate of ionizing 
photons (his > 13.6eV) emitted by the secondary star per steradian, N^- F° r the secondary 
of 7] Car N i2 = $ 2 /47r = 2 x 10 48 s~ 1 sr~ 1 , as derived by Kashi & Soker (2007a) based on 
Schaerer & de Koter (1997). We neglect recombination in the low density secondary wind. 
Thus the distance along the line from the contact discontinuity through the shocked primary 
wind where the wind is optically thin to the secondary's radiation, as can be seen in Fig. 
EH is 

<4bs = 9 (Bl) 

where is the recombination coefficient, r 2 is the distance to the secondary, n e is the 
electron number density and uh is the proton number density. 



This preprint was prepared with the AAS IATgX macros v5.2. 
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Fig. 9. — The three stages of the conical shell ionization plotted along the secondary's orbit 
together with the orbital phases of transitions between stages. Darker colors means denser 
wind. Stage (i): When the orbital separation is large, the hot secondary ionizing radiation 
manages to keep the conical shell (the post-shock primary wind zone) totally ionized. In 
that stage dust can be formed only at very large distances. Stage (ii): When the orbital 
separation decreases to r < 14.3 AU at ~ 170 days (9 = 150°) before (after) periastron, 
the density of the gas in the conical shell increases, the cone becomes optically thick to the 
secondary's UV radiation and parts of the cone begin to be neutral. First at large distances, 
and then closer and closer to the stagnation point along the line joining the two stars. Stage 
(Hi): As the secondary approaches periastron, r < 3.2 AU at ~ 20 days (9 = 90°) before 
(after) periastron the ionization front is within the conical shell in all directions, and the 
region near the stagnation point becomes neutral. During most of this stage the conical 
shell does not exist according to our model. Instead, an accretion flow onto the secondary 
takes over. The exact times (and locations) of the transitions between the stages depend on 
poorly determined winds' parameters. 

We recall that the ionizing radiation from the primary is neglected here because most 
of it is absorbed in the dense gas of the primary wind close to the primary, and in the dense 
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gas in the conical shock region. Despite being five times as luminous as the secondary, the 
primary is much cooler and its ionizing photon flux (hv > 13.6 eV) is small to begin with. 
However, at longer wavelengths in the UV and in the visible band the primary radiation 
cannot be neglected, and it is considered as a heating source of the dust together with the 
secondary radiation (see discussion following equation [7]). 

The value of d a b s (a) is then compared to the thickness d(a) of the wind in the same 
direction a. In stage (i) <i a bs > d for every a and the shocked primary wind (the conical 
shell) is totally ionized. In stage (ii) d^ > d for small values of a and d a ^ s < d for large 
values of a, and parts of the cone begin to be neutral. In stage (Hi) d a \> s < d for every a 
and the ionization front is within the conical shell in all directions. These three stages are 
drawn schematically in Figs. 151 and [HH 

When the material in the conical shell is neutral, the conditions are favorable for dust 
to be formed. The gas in the conical shell is initially hot T gas = 10 4 K. As it expands and 
cools, hot dust forms at temperatures T < 1700 K (Van Genderen et al. 1994). The cooling 
time depends on density. We use figure 11 of Woitke et al. (1996) to estimate the cooling 
time from T = 10 4 K and n ~ 10 12 — 10 14 cm" 3 to T = 1700 K in an adiabatic process and 
find it to be approximately 5 — 10 days. 

For most of the orbital period the wind is ionized and no dust is being formed in the 
conical shell close to the secondary. The roughly constant IR emission is emitted by dust 
which was formed in previous cycles, mostly cooler dust, and the hotter dust that is formed 
in the primary wind and heated by the secondary and primary radiation. The rapid increase 
in the IR emission ~ 150 days before the 2003.5 event (see section [7]) can be explained by 
the process described in stage (ii) - more and more material in the conical shell becomes 
neutral, favoring the formation of dust closer to the secondary. The occurrences of the event 
are explained in section [7J 

We stress again that the numbers can be somewhat different for other parameters, like 
e, stellar masses and magnetic field in the primary wind. 
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Fig. 10. — Schematic drawing of the plane between the two stars, at the three stages of 
the conical shell ionization, d is the thickness of the shocked primary wind zone (conical 
shell) from the contact discontinuity to the primary shockwave in a direction. <i a bs is the 
distance along the line from the contact discontinuity through the shocked primary wind 
where the wind is optically thin to the secondary's radiation, in the same direction. Upper 
panel — Stage(i): Until ~ 170 days before (after) periastron <i abs > d for all values of a. The 
secondary ionizing radiation manages to keep the conical shell totally ionized. Middle panel 
— Stage(ii): From ~ 170 days until ~ 20 days before (after) periastron rf abs > d for small 
values of a and <i a bs < d for large values of a, and parts of the cone begin to be neutral. 
Lower panel — Stage(m): From ~ 20 days before periastron until ~ 20 days after periastron 
^abs < d for every a and the ionization front is within the conical shell in all directions, and 
the region near the stagnation point becomes neutral. The last stage exists only for a short 
period of time, or might not exist at all, because of the collapse of the conical shell onto the 
secondary ~ 20 days before periastron, and the formation of an accretion flow (Akashi et al. 
2006). 
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Fig. 11. — dab s (a) is the distance along the line at angle a from the contact discontinuity 
through the shocked primary wind (conical shell) required to absorb the ionizing radiation 
of the secondary star. In our model the asymptotic angle of the conical shell is at a = 60° 
(left side of the figure), where a = 180° is the angle toward the primary and a = 0° in a 
direction opposite to the primary. d(a) is the total length through the conical shell along the 
same direction (see Fig. [TP]) . If d^^a) / d(a) > 1 then the conical shell along the direction a 
is completely ionized, while if d^^a) / d{a) < 1 part of the conical shell is neutral. Neutral 
regions are more favorable for dust formation. 9 is the orbital angle (see Fig. [9]), with 6 = 
at periastron (t = 0), and 6 = 180° at apastron (t = 1012 days). The three stages discussed 
in the text, and the times (in days) before periastron passage of transitions between the 
stages are marked on the figure. 



